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Dynamic  Stability  of  Maglev  Systems 

Y.  Cai,  S.  S.  Chen,  T.  M.  Mulcahy,  and  D.  M.  Rote 


Abstract 


Because  d3mamic  instability  is  not  acceptable  for  any  commercial  maglev 
systems,  it  is  important  to  consider  this  phenomenon  in  the  development  of  all 
maglev  systems.  This  study  considers  the  stability  of  maglev  systems  based  on 
experimental  data,  scoping  calculations,  and  simple  mathematical  models. 
Divergence  and  flutter  are  obtained  for  coupled  vibration  of  a  three-degree-of- 
freedom  maglev  vehicle  on  a  guideway  consisting  of  double  L-shaped  aluminum 
segments  attached  to  a  rotating  wheel.  The  theory  and  analysis  developed  in  this 
study  identifies  basic  stability  characteristics  and  future  research  needs  of  maglev 
systems. 


1  Introduction 


The  dynamic  response  of  maglev  systems  is  important  in  several  respects: 
safety  and  ride  quality,  guideway  design,  and  system  costs.  Ride  quality  is 
determined  by  vehicle  response  and  by  environmental  factors  such  as  humidity 
and  noise.  The  dynamic  response  of  vehicles  is  the  key  element  in  the 
determination  of  ride  quality,  and  vehicle  stability  is  one  of  the  important 
elements  relative  to  safety.  To  design  a  proper  guideway  that  provides  acceptable 
ride  quality  in  the  stable  region,  the  vehicle  dynamics  must  be  understood. 
Furthermore,  the  trade-ofT  between  guideway  smoothness  and  the  levitation  and 
control  systems  must  be  considered  if  maglev  systems  are  to  be  economically 
feasible.  The  link  between  the  guideway  and  other  maglev  components  is  vehicle 
d3mamics.  For  a  commercial  maglev  system,  vehicle  dynamics  must  be  analyzed 
and  tested  in  detail. 

For  safety,  maglev  systems  should  be  stable.  Magnetic  forces  are  basically 
position-dependent,  although  some  are  also  velocity-dependent.  These  motion- 
dependent  magnetic  forces  can  induce  various  t3rpes  of  instability.  In  addition, 
the  periodic  structure  of  the  motion-dependent  magnetic  forces  may  in  some  cases 
also  induce  parametric  and  combination  resonances. 
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Some  analytical  and  experimental  studies  have  been  performed  to 
understand  the  stability  characteristics  of  maglev  systems.  Several  examples  are 
summarized  briefly  as  follows: 

1.1  Theoretical  Studies 

*  Davis  and  Wilkie  (1971)  studied  a  magnetic  coil  moving  over  a 
conducting  track  and  concluded  that  negative  damping  occturs  at 
velocities  greater  than  the  characteristic  velocity  based  on  thin-track 
theory. 

*  Ohno  et  al.  (1973)  studied  the  pulsating  lift  forces  in  a  linear 
synchronous  motor.  These  forces  may  cause  parametric  and 
combination  resonances,  as  well  as  heave  and  pitch  oscillations. 

*  Baiko  et  al.  (cited  in  Chu  and  Moon  1983)  considered  the  interactions 
of  induced  eddy  currents  with  on-board  superconducting  magnets 
and  found  possible  heave  instabilities. 

1.2  Exparlmental  Studies 

*  An  experimental  vehicle  floating  above  a  large  rotating  wheel  was 
found  by  Moon  (1974)  to  have  sway-yaw  instabilities. 

*  Experiments  performed  at  MIT  on  a  test  track  showed  pitch-heave 
instability  (Moon  1975). 

1.3  Experlmental/Analytical  Studies 

*  A  conducting  guideway,  consisting  of  L-shaped  aluminum  segments 
attached  to  a  rotating  wheel  to  simulate  the  Japanese  full-scale 
guide  way  at  Miyazaki,  was  studied  experimentally  and  analytically 
by  Chu  and  Moon  (1983).  Divergence  and  flutter  were  obtained  for 
coupled  yaw-lateral  vibration;  the  divergence  leads  to  two  stable 
equilibrium  yaw  positions,  and  the  flutter  instability  leads  to  a  limit 
cycle  of  coupled  yaw  and  lateral  motions  in  the  neighborhood  of  the 
drag  peak. 

*  Variation  of  the  magnetic  lift  force  due  to  variation  of  the  levitated 
height  corresponding  to  the  sinusoidal  guideway  roughness  was 
studied  by  Yabuno  et  al.  (1989).  Parametric  resonance  of  heaving  and 
pitching  motions  is  possible. 
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Based  on  these  published  analytical  results  and  experimental  data,  it  is 
obvious  that  different  types  of  d3mamic  instabilities  can  occur  in  maglev  systems. 
Because  dynamic  instability  is  not  acceptable  for  any  commercial  maglev 
systems,  it  is  important  to  consider  this  phenomenon  in  the  development  of  all 
maglev  systems. 

This  study  considers  the  stability  of  maglev  systems  and  is  based  on 
experimental  data,  scoping  calculations,  and  simple  mathematical  models.  The 
objective  is  to  provide  some  basic  stability  characteristics  and  to  identify  future 
research  needs. 


2  Motion-Dependent  Magnetic  Forces _ 

2.1  Motion-Dependent  Magnetic-Force  Coefficients 

Magnetic  forces  are  needed  for  any  vehicle  dynamics  analysis,  guideway 
structural  design,  fastening  design,  and  prediction  of  ride  quality.  These  force 
components  are  considered  from  the  standpoint  of  vehicle  stability. 

As  an  example,  consider  a  vehicle  with  six  degrees  of  freedom,  three 
translations,  Ux,  Uy,  Uj,  and  three  rotations,  ©x,  ©y,  ©z,  as  shown  in  Pig.  1.  Let  U 
be  the  vector  consisting  of  the  six  motion  components,  i.e.. 
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Velocity  and  acceleration  are  given  by 


and 


U  = 


(2) 


The  motion-dependent  magnetic  forces  can  be  written 
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Fig.  1.  Displacement  components  of  a 
maglev  system 
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j=l 

where  mij,  qj,  and  ky  are  magnetic  mass,  damping,  and  stiffness  coefficients. 
These  coefficients  can  be  obtained  anal3rtically,  numerically,  or  experimentally 
and  are  functions  of  the  system  parameters. 

*  Analytical  Studies:  Analyses  for  simple  cases  can  be  performed  to 
determine  the  characteristics  of  the  coefficients.  For  example,  an 
analytical  method  may  be  used  to  identify  the  coefficients  that  can  be 
neglected  under  specific  conditions. 

*  Numerical  Methods:  For  the  general  case  with  complicated 
geometries,  analytical  methods  may  not  be  appropriate  and 
numerical  methods  will  be  more  useful.  Numerical  methods  (finite- 
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element  method  and  boundary-element  method)  can  be  used  to 
calculate  the  values  of  all  coefficients  under  specific  conditions. 

*  Experimental  Techniques:  Measurements  of  magnetic  forces  will 
give  the  information  required  to  calculate  magnetic-force  coefficients. 


2.2  Experimental  Methods  to  Measure  Motion-Dependent 
Magnetic-Force  Coefficients 

Quasistatic  Motion  Theory.  The  magnetic  forces  acting  on  an  oscillating 
vehicle  are  equal,  at  any  instant  in  time,  to  those  of  the  same  vehicle  moving  with 
a  constant  velocity  and  with  specific  clearances  equal  to  the  actual  instantaneous 
values.  The  magnetic  "orces  depend  on  the  deviation  from  a  reference  state  of 
speed  and  clearance,  i.e.,  the  motion-dependent  magnetic  forces  depend  only  on 
Uj,  but  not  Uj  and  iij,  so  that 

6 

U  =  ZkyUj.  (4) 

j=l 


In  this  case,  the  magnetic  forces  are  determined  uniquely  by  vehicle  position.  All 
elements  of  magnetic  stiffness  kjj  can  be  obtained.  To  determine  kjj,  the 
magnetic-force  component  fj  is  measured  as  a  function  of  uj.  Stiffiiess  is  given  by 


duj 


(5) 


In  general,  kjj  is  a  function  of  U. 

Unsteady-Motion  Theory.  The  magnetic  forces  acting  on  an  oscillating 
vehicle  will  depend  on  U,  U,  and  U.  The  magnetic  force  based  on  the  unsteady- 
motion  theory  can  be  obtained  by  measuring  the  magnetic  force  acting  on  the 
vehicle  oscillating  in  the  magnetic  field.  For  example,  if  the  displacement 
component  uj  is  excited,  its  displacement  is  given  by 

u  j  =  u  j  exp(V^  tot).  (6) 

The  motion-dependent  magnetic  force  of  the  component  fj  acting  on  the  vehicle  is 
given  by 
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fj  =  cos(\|fy )  +  ay  sin(vy  )ju  j  exp^V^  cot),  (7) 

where  ay  is  the  magnetic  force  amplitude  and  is  the  phase  angle  between  the 
magnetic  force  and  the  vehicle  displacement  uj.  These  values  are  measured 
experimentally. 

Using  Eqs.  3  and  6,  we  can  also  write  the  motion-dependent  magnetic  force 
component  as 

fj  =  ^-niyto^  +  (0  cy  +  )u  j  exp(V^  cot).  (8) 

A  comparison  of  Eqs.  7  and  8  yields 

Cij  =  aijSin(\|/y)/co, 

(9) 

my=[ky-aycos(vy)]/co2. 

Based  on  Eqs.  5  and  9,  all  motion-dependent  magnetic-force  matrices  can  be 
determined  from  two  experiments:  quasistatic  motion  and  imsteady  motion. 

If  mij  and  cy  are  of  no  concern,  the  experiment  using  quasisteady  motion  is 
sufficient  to  determine  ky. 


2.3  Quasistatic  Motion-Dependent  Magnetic-Force  Coefficients  of 
Magiev  System  with  L-Shaped  Guideway 

An  experiment  w  s  conducted  recently  at  Argonne  National  Laboratory  to 
investigate  the  lift,  drag,  and  guidance  magnetic  forces  on  an  NdFeB  permanent 
magnet  moving  over  an  aluminum  (6061-T6)  L-shaped  ring  mounted  on  the  top 
surface  of  a  1.2-m  diameter  rotating  wheel  (.shown  in  Fig.  2).  For  a  given  rotating 
speed  of  the  wheel,  the  lift  and  guidance  magnetic  forces  were  measured  as  the 
guidance  gap  Y*  and  lift  height  h  were  varied.  A  schematic  diagram  of  the 
measurement  approach  is  shown  in  Fig.  3.  Figure  4  shows  those  measured 
forces  as  a  function  of  h,  with  Y*  fixed  (Y*  =  5  mm  and  12.7  mm),  or  as  a  function 
of  Y*  with  h  fixed  (h  =  7  mm  and  12.7  mm)  when  the  surface  velocity  of  the  lateral 
leg  of  the  guideway  is  36.1  m/s,  the  highest  velocity  tested. 

During  testing,  the  long  side  of  the  25.4  x  50.8  x  6.35  mm  rectangular  magnet 
was  oriented  parallel  to  the  direction  of  motion  of  the  L-shaped  guideway  and  was 


Fig.  2.  Experimental  apparatus  for  magnetic  force  measurement 


Fig.  3.  Schematic  diagram  of  the  apparatus  used  to  measure 
magnetic  forces  on  an  L-shaped  aluminum  sheet 
guideway 
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held  stationary  by  a  two-component  force  transducer  that  comprised  two  BLH 
C2G1  load  cells  connected  in  series  to  measure  the  lift  and  guidance  forces 
simultaneously.  Laboratory  weights  were  used  to  calibrate  the  transducer  and  to 
assess  crosstalk  (which  was  found  to  be  <2%).  The  base  of  the  load  cell  assembly 
was  mounted  on  motorized  stages  that  provided  accurate  positioning  (±0.05  mm). 
Out-of-roundness  of  the  L-shaped  gmdeway  ring  varied,  but  was  always  less  than 
±0.15  mm  for  the  lateral  leg  and  ±0.35  mm  for  the  vertical  leg.  Ability  to  exactly 
position  the  magnet  with  respect  to  the  guideway  dominated  our  experimental 
error,  estimated  at  ±5%. 

The  qualitative  trends  in  the  lift  and  guidance  force  data  taken  at  lower 
velocities  for  the  L-shaped  guideway  are  the  same  as  shown  in  Fig.  4  for  the 
highest  velocity  tested.  Two  interesting  features  are  evident.  First,  a  maximum 
occurs  in  the  guidance  force  variation  with  respect  to  height  variations  at  a  fixed 
gap  (as  shown  in  Figs.  4e  and  4f)  that  is  caused  by  the  comer  region  of  the 
guideway.  Second,  a  minimum  occurs  in  the  guidance  force  variation  with 
respect  to  gap  variations  at  a  fixed  height  (as  shown  in  Figs.  4g  and  4h)  that  is 
caused  by  the  edge  of  the  lateral  leg  of  the  guideway.  As  will  be  shown  in  Section 
4.4,  the  first  feature  is  associated  with  a  flutter  and  the  second  with  a  divergence 
instability. 

Based  on  the  magnetic  force  data  shown  in  Fig.  4,  we  can  calculate  the 
quasistatic  motion-dependent  magnetic-force  coefficients  with  Eq.  5.  All  elements 
of  magnetic  stiffness  k//,  k/g,  kg/,  and  kgg,  were  calculated  and  are  shown  in 
Fig.  5  with  various  Y*  and  h. 

The  curve  fits  to  both  magnetic  forces  and  stiffnesses  were  derived  with 
pol3momial  expressions  (results  are  given  in  Figs.  4  and  5)  and  input  into  a 
computer  code  to  simulate  coupled  vibrations  of  the  maglev  vehicle. 


3  Stability  of  Maqiev  Systems _ 

Without  motion-dependent  magnetic  forces,  the  equation  of  motion  for  the 
vehicle  can  be  written  as 

[M,]{0}+(C,]{u}+[K,]{U}  =  {Q},  (10) 

where  My  is  the  vehicle  mass  matrix,  Cy  is  the  vehicle  damping  matrix,  Ky  is  the 
vehicle  stiffness  matrix,  and  Q  is  the  generalized  excitation  force. 

The  motion-dependent  magnetic  forces  are  given  in  Eq.  3.  With  motion- 
dependent  magnetic  forces,  Eq.  10  becomes 
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[M,  +  M„]{U}  +  [C,  +  C„]{U}  +  [K,  +  K„]{U)  =  {Q}.  (11) 

where  Mjq  is  the  magnetic  mass  matrix,  is  the  magnetic  damping  matrix, 
and  Kjq  is  magnetic  stiffness;  their  elements  are  m^,  Cij,  and  k^. 

Once  the  magnetic-force  coefficients  are  known,  analysis  of  vehicle  stability 
is  straightforward.  Equation  11  may  be  written  as 

[M]{U}+[C]{U}+[K]{U}  =  {Q}.  (12) 

In  general,  M,  C,  and  K  are  functions  of  U,  U,  and  0;  therefore,  a  complete 
solution  is  rather  difficult  to  obtain.  In  many  practical  situations,  one  can  ignore 
all  nonlinear  terms,  so  that  M,  C,  and  K  are  independent  of  vehicle  motion. 

By  premultiplying  by  {U}'*’  and  forming  the  S3rmmetric  and  antisymmetric 
components  of  the  matrices 

[Ml]  =  i((M) +[M]’'),  [Mj]  =  |([M]  -  (Mf), 

[Ci]  =  i(lC]+(Cf).  [C2]  =  |([C)-tC]’'),  (13) 

[KJ  =  |([K]+ [Kf ),  [K2I  =  |(tK)  -(Kf ), 

the  terms  may  be  separated,  giving 

{uf[Mi]{u}  +  {ijf  [C2]{U}+ {Uf  [KiKU] 

=  -[{l7}’'[M2]{u}  +  {uf[Ci]{u}+{u}’'[K2KU})  (14) 

Equation  14  equates  rates  of  work.  The  terms  on  the  right-hand  side  of  the 
equation  produce  a  net  work-resultant  when  integrated  over  a  closed  path 
through  the  space  {U},  the  magnitude  depending  on  the  path  taken.  The  forces 
corresponding  to  the  matrices  M22  Ci,  and  K2,  appearing  on  the  right-hand  side, 
are  thus,  by  definition,  the  nonconservative  parts  of  the  forces  represented  by  M, 
C,  and  K.  The  terms  on  the  left-hand  side  similarly  can  be  shown  to  give  rise  to  a 
zero  work-resultant  over  any  closed  path,  and  therefore  together  are  the  sum  of 
the  rates  of  work  from  the  potential  forces  and  the  rate  of  change  of  kinetic  energy. 
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Different  tjrpes  of  instability  can  be  classified  according  to  the  dominant  terms  in 
Eq.  14: 

•  Magnetic-Damping'Controlled  Instability  (single-mode  flutter):  The 
dominant  terms  are  associated  with  the  symmetric  damping  matrix 
[Cil.  Flutter  arises  because  the  magnetic  damping  forces  create 
"negative  damping,"  that  is,  a  magnetic  force  that  acts  in  phase  with 
the  vehicle  velocity. 

•  Magnetic-Stiffness-Controlled  Instability  (coupled-mode  flutter):  The 
dominant  terms  are  associated  with  the  antisymmetric  stiffness 
matrix  [K2].  It  is  called  coupled-mode  flutter  because  at  least  two 
modes  are  required  to  produce  it. 

Corresponding  to  the  single-  and  coupled-mode  flutter,  parametric  and 
combination  resonances  may  exist  if  the  motion-dependent  magnetic  forces  are  a 
periodic  function  of  time. 

•  Parametric  Resonance:  When  the  period  of  a  motion-dependent 
magnetic  force  is  a  multiple  of  one  of  the  natural  frequencies  of  the 
vehicle,  the  vehicle  may  be  d3mamically  unstable. 

•  Combination  Resonance:  When  the  period  of  the  motion-dependent 
magnetic  forces  is  equal  to  the  sum  or  difference  divided  by  an 
integer  of  the  natural  frequencies  of  the  vehicle,  the  vehicle  may  also 
be  subjected  to  dynamic  instability. 

In  practical  cases,  two  or  more  mechanisms  may  interact  with  one  another,  and 
Eq.  12  is  applicable  for  general  cases. 

It  is  noted  that  maglev  systems  are  subjected  to  several  groups  of  forces, 
including  magnetic  forces,  aerodynamic  forces,  and  guideway  perturbation.  The 
theory  presented  in  this  paper  is  applicable  to  maglev  systems  when  they  are 
subjected  to  other  tjrpes  of  forces.  In  partietdar,  the  aerodynamic  effects  can  be 
described  exactly  the  same  way  as  those  given  in  Eqs.  1-14  and  the  dynamic 
response  characteristics  to  aerodynamic  forces  are  similar  to  magnetic  forces;  see 
Chen  (1987)  for  details. 


4  Simplified  Vehicle  Models  for  Dynamic  Instability _ 

Four  different  vehicles  are  considered,  in  order  to  provide  an  understanding 
of  stability  characteristics. 
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4..  Two-Degree>of-Freedom  Vehicle 

A  maglev  vehicle  is  supported  by  magnetic  forces;  the  resultant  lift  and  drag 
forces  of  the  coil  above  the  continuous*sheet  track  can  be  represented 
approximately  by  (Sinha  1987) 

FL(v,z,t)  =  eLF(t), 

FD(v,z,t)  =  SoFft), 

F(t)  =  PoI2/47cz, 

(15) 

El  =  1  -  1/(1  +  v^/w^)”, 

Ejj  =  (w/v)[l  -  1/(1  +  v2/w2)"], 
w  =  2/poOh, 

where  v  is  the  forward  velocity,  z  is  the  steady-state  height  of  the  coil  above  the 
track,  t  is  the  time,  I  is  the  constant  coil  cxirrent,  and  w  is  the  characteristic  speed 
and  is  related  to  track  thickness  h,  conductivity  o,  and  permeability  Po.  The  value 
of  n  is  1  for  a  single  conductor  and  varies  from  about  1/5  to  1/3  for  coils  (Rhodes 
and  Mulhall  1981).  Note  that  F(t)  represents  the  repulsive  force  between  the  coil 
and  its  image  coil.  The  force  ratios,  El  andep,  are  given  in  Fig.  6  for  n  s  1,  1/3, 
and  1/5. 


Assume  that  the  vehicle  is  traveling  at  a  velocity  Vg  at  cm  equilibrium  height 
Zq.  The  instantaneous  position  and  height  of  the  vehicle  are  x  and  z  respectively; 
therefore, 

x(t)  =  Vot  +  X(t), 

(16) 

Z(t)  =  Zq  +  Z(t). 

The  equations  of  motion  for  the  vehicle  moving  at  a  velocity  v(t)  with  a  levitation 
height  z(t)  can  be  written  as 


m2(t)  a  -mg  +  FL(v,z,t), 
m  jt(t)  s  Fp  -  FD(v,z,t)  -  Fa, 


(17) 


where  m  is  the  mass  of  the  vehicle,  Fp  is  the  propulsion  force,  and  Fa  is  the  aero¬ 
dynamic  force.  The  aerodynamic  force  is  given  by 


20 


1 


0.8 

0 


4 

0.2 


0 

0  10  20  30  40  SO 

V/W 

Fig.  6.  Magnetic  forces  divided  by  image  force 

F,  =  K,y2 

s  0.5  Cj)LAp, 


(18) 


where  L  and  A  are  the  length  and  cross-sectional  area  of  the  vehicle,  p  is  air 
density,  Cd  is  the  drag  coefficient  (with  its  values  varying  from  0.2  to  0.3, 
depending  on  vehicle  shape). 


The  equilibrium  point  of  the  vehicle  v^  and  Zq  is  defined  by 


Fl(Vo,2o)  =  “g. 
FD(Vo,Zo)  =  Fp-Fa- 


(19) 


Using  Eqs.  16-19  and  neglecting  the  nonlinear  terms,  we  obtain  the  following 
equations  of  motion  for  the  vehicle  X(t)  and  Z(t): 


m2(t)  -  C2x(v,z,t)X(t)  +  Kx^(v,z,t)Z(t)  =  0, 
m  jt(t)  +  Cxx(v,z,t)X(t)  -  Kjtj5(v,z,t)Z(t)  =  0 


(20) 
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CjxCvjZ.t)  and  Cxx(v,z,t)  are  magnetic  damping,  and  Kz2(v,z,t)  and  Kxz(v,z,t)  are 
magnetic  stiffness;  these  are  given  by 

Cxx(v,z,t)  =  aFL(v,z,t)/av  I  (v^,,Zo)» 

Cxx(v,z,t)  =  aFD(v,z,t)/9v  I  (v^„Zo)  +  2KaVo, 

Wl/ 

Kzz(v,z,t)  =  -  aFL(v,z,t)/3z  I  (Vo,Zo)» 

Kxz(v,z,t)  =  -aFD(v,z,t)/az  I  (Vo,Zo)* 

Their  values  are  given  as  follows; 

Czx  =  (mg/eL)2nvo/[w2(l  +  vj/w2)n+l], 

Cxx  =  wmg/vj  [w2  -  (2n  -  I)v2/(w2  +  vj)]  +  2KaVo, 

(22) 

Kzz  =  mg/zo, 

Kxz  =  wmg/voZQ. 

For  high  speed  vehicles,  the  values  of  Czx  and  Kxz  approximately  zero. 
Therefore,  the  motions  in  the  vertical  direction  and  forward  direction  are 
imcoupled  at  high  speeds. 

Equations  20  can  be  anedyzed;  let 

Z(t)  =  a  exp(i(it)t), 

(23) 

X(t)  =  b  exp(io)t). 

Substituting  Eqs.  23  into  Eqs.  20  gives  the  following  frequency  equation: 


-mcD^  +  Kjjx  -iojCxx 


a1  fO 


-mco^  +  icoCxxJ  lOJ 


The  natural  frequencies  can  be  determined  from  the  determinant  of  the  coefficient 
matrix  given  in  Eq.  24.  At  high  speeds,  the  off-diagonal  terms  may  be  neglected. 
The  natural  frequency  of  the  vertical  motion  fy  is 

fy  =  (g/Zo)0-5/2jt.  (25) 
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Fig,  7.  Natural  frequency  as  a  function  of  levitation  height 

The  natural  frequency  in  the  vertical  direction  depends  on  the  levitation  height 
only.  Figure  7  shows  the  natural  frequent^  as  a  function  of  the  levitation  height. 
At  high  speeds,  oscillations  in  the  vertical  direction  are  stable. 

In  the  forward  direction,  the  motion  is  given  by 

X(t)  =  Cl  +  C2exp(st).  (26) 

For  high-speed  vehicles,  the  exponent  s  is  approximately  given  by 

s  =  (2n  -  l)wg/v2  -  2KaVo.  (27) 

Note  that  the  vehicle  may  be  unstable  if  n  =  1,  and  is  zero.  At  high  speeds,  the 
second  term  given  in  Eq.  27  is  larger  than  the  first  term  regardless  of  the  values  of 
n;  therefore,  s  is  negative  and  the  system  is  stable. 


4.2  Three-Degree-of-Freedom  Vehicle 

Figure  8  shows  a  three-degree-of-freedom  vehicle  traveling  at  a  velodiy  Vo  at 
an  equilibrium  height  Zq.  For  a  S3rmmetric  vehicle,  the  instantaneous  position 
and  height  of  the  vehicle  are  x(t),  zi(t),  and  Z2(t);  therefore. 
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Fig.  8.  Three-degree-of-freedom  vehicle 
x(t)  =  Vo(t)  +  X(t) , 

zi(t)  =  Zo  +  Zi(t) ,  (28) 

Z2(t)  =  Zo  +  Z2(t) . 

The  equations  of  motion  for  the  vehicle  moving  at  a  velocity  v(t)  with  levitation 
height  zi(t)  and  Z2(t)  can  be  written  as 

mx(t)  =  Fp  -FDi(zi,v,t)-FD2(z2,v,t)-Fa, 

^[zi(t)  +  Z2(t)]  =  -mg  +  FLi(zi,v,t)  +  Fl2(z2,  v,t),  (29) 

^[zi(t)  -  Z2(t)]  =  ^[FLi(zi,v,t)  -  Fl2(z2,  v,t)]. 

where  m  is  the  vehicle  mass,  le  is  the  rotational  moment  of  inertia  about  the 
vehicle's  center  of  mass,  Fp  is  the  propulsion  force,  and  Fa  is  the  aerodynamic 
force  that  is  assumed  to  act  at  the  center  of  mass  of  the  vehicle  and  is  given  by 
Eq.  18. 
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For  a  symmetric  vehicle  with  two  identical  levitation  systems  at  the  two  ends, 
the  eqmlibrium  point  of  the  vehicle  Vq,  ziq,  and  Z20>  as  well  as  the  magnetic  forces, 
is  defined  as 

Fli(Vo,zio)  +  Fl2(Vo,Z20)  =  mg  , 

Pdi(Vo,zio)  +  Fd2(Vo,Z20)  =  Fp  -  Fa , 

zio  =  Z20  =  zo ,  (30) 

Fli(vo,zio)  =  Fl2(vo,Z20)  , 

Fdi(Vo,Z20)  ”  Fd2(Vo,Z20)  • 

Using  Eqs.  28,  29,  30,  and  18  and  neglecting  the  nonlinear  terms,  we  obtain  the 
following  equations  of  motion  of  the  vehicle,  X(t),  Zi(t),  and  Z2(t): 

Y(Zi  +  Z2) + (Cjxi  +  Czx2)X  +  K„iZi  +  K222Z2  =  0 

^(zi  -  Z2) + (Czxi  -  Czx2)X  +  (K„iZi  -  K„2Z2  =  0  (31) 

mX  +  (Cjjx  1 + Cxx2  +  ^kgVo  )X  +  Kx2iZi  +  K„2Z2  -  (^ 

where 


i  =  1, 2. 
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The  magnetic  damping  coefficients  Czxi  £^d  Cxxi>  and  the  magnetic  stiffness 
coefficients  Kzzi  and  Kxzi  can  he  calculated  from  the  magnetic  lift  and  drag  forces 
given  in  Eq.  15. 

At  high  speeds,  Kzxi  and  Cxzi  are  approximately  zero.  The  equations  of 
motion  become 


|(Z,  +  Z2)+K,,iZi  +  K,,2Z2=0, 

^(Zi-Z2)  +  K„iZi-K,,2Z2  =  0.  (33) 

Li 


mX  +  (Cxxl  +  Cxx2  +  2KaVo)X  =  0. 


In  this  case,  the  vehicle  is  stable  at  high  speeds;  this  is  similar  to  the  two-degree- 
of-freedom  vehicle.  The  natural  frequency  of  vertical  oscillations  is  the  same  as 
that  in  the  two-degree-of-freedom  system  given  in  Eq.  25.  The  natural  frequency  of 
pitching  oscillations  is 


27C 


4l02o 


(34) 


For  a  square  vehicle  with  length  L  and  height  h,  the  natural  frequency  of 
pitching  oscillations  is 


2n 


/  > 

0.5 

f  \ 

g 

3 

i+% 

>.0.5 


(35) 


The  natural  frequency  of  pitching  oscillations  is  larger  than  vertical  oscillations. 
For  a  long  vehicle  (h  «  L),  fp  is  equal  to  about  1.7  fy.  For  a  square  vehicle,  fp  =  1.4 
fy.  At  high  speeds,  heaving  and  pitching  oscillations  are  stable  for  the  magnetic 
levitation  described  by  Eq.  15.  In  the  forward  direction,  the  result  is  the  same  as 
that  for  the  two-degree-of-freedom  vehicle. 
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4.3  Six-Degree-of-Freedom  Vehicle 

For  the  six-degree-of-freedom  vehicle  shown  in  Fig.  1,  stability  can  be  studied 
from  Eqs.  11  or  12.  Once  the  coefficients  of  magnetic  forces  m^,  Cij,  and  ky  are 
known,  Eqs.  11  or  12  can  be  evaluated  for  Q  =  0.  Let  the  displacement  of  a 
particular  component  be 

uj(t)  =  ajexpCX  +  icoH  .  (36) 

Substituting  Eq.  36  into  Eq.  12  with  Q  =  0  yields 

{(X  +  iG))2[M]  +  (X  +  i(o)[C]  +  [K]}  {A}  =  {0} .  (37) 

The  values  of  X  and  o)  can  be  calculated  based  on  Eq.  37  by  setting  the  determinant 
of  the  coefficient  matrix  equal  to  zero. 

Vehicle  stability  is  determined  by  X,  which  is  a  function  of  v.  If  X  <  0,  vehicle 
motion  is  damped;  if  X  >  0,  vehicle  displacement  increases  with  time  until 
nonlinear  effects  become  important. 

To  solve  this  problem,  all  motion-dependent  magnetic-force  coefficients  must 
be  known.  At  this  time,  it  appears  that  limited  analytical,  numerical,  or 
experimental  data  are  available.  For  any  future  maglev  systems,  it  will  be 
necessary  to  investigate  the  characteristics  of  motion-dependent  magnetic  forces 
to  avoid  dynamic  instabilities. 


4.4  Vehicle  on  Double  L-Shaped  Aluminum  Sheet  Guideway 

Figure  9  shows  a  cross  section  of  a  vehicle  on  a  double  L-shaped  aluminum 
sheet  guideway.  Assume  that  the  vehicle  travels  at  a  constant  velocity  along  x 
direction.  Two  permanent  magnets  are  attached  to  the  bottom  of  vehicle  and 
provide  lift  and  guidance  force  Fl^,  Fl2>  Fgj,  and  Fg2  (see  Fig.  9).  Assume  at  the 
initial  state  that  hj  s  h2  -  ho  and  gi  =  g2  =  go>  thus,  the  vehicle  and  guideway 
geometries  can  be  expressed  as  follows: 

Li  s  L2  -  S  s  76.2  (mm) 

W  =  152.4  +  S  -  2go  (mm) 


H=:0.9W  (mm) 
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Fig.  9.  Maglev  system  with  a  vehicle  operating  on  double  L- 
shaped  aluminum  sheet  guideway 

a  s  0.5  H  (mm) 

b  ss  0.5(W  -  25.4)  (mm). 

Equations  of  motion  for  this  three-degree-of-freedom  maglev  system  can  be 
written  as 

mz  +  Cz  =  FLi+FL2  -mg 


10 + E0  =  (Fqj  +  Fg2  )a  +  (Fqj  +  Fg2  )b 


(38) 


my  +  Dy  =  FGi+FG2, 


where  m  is  the  mass  of  the  vehicle,  C  and  D  are  damping  ratios;  I  is  the  moment 
of  inertia  about  the  center  of  mass  inertia  moment  of  the  vehicle  [I  s 
(m/12)(H2  +  W2)].  Flj,  Fl2,  Fgj,  and  Fg2  are  lift  and  guidance  forces  and  are 
functions  of  y  and  z.  At  the  equilibrium  position,  they  are  FLio(yo>Zo)>  Fl2o^o>Zo)> 
FGio(yo.Zo)»  and  FG2o(yo,Zo)-  Apply  them  to  Eqs.  38: 


Flio*=FL0o 
FLio  +  FL20  =  mg 
FQio  -  “FG20 
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(39) 


Therefore 

^  2F,.(h„) 

g  g 


(40) 


Let 


z 


|(«l  +  «2) 


y=u3 


0  =  (ui-U2)/2b 


(41) 


where  ui,  U2,  and  U3  are  shown  in  Fig.  10.  Equations  38  can  be  rewritten  as 
m(tii  +  02) + c(ui  -  U2)  =  2(Fli  +  Flj  -  “g) 

^(ui-U2)+^(ui-U2)  =  2a(FGi  +FG2)+2b(FGj  -FGg)  (42) 

mii3  +  Du3  =  Fgj  +  Fgj  • 

Note  the  reduced  dependence  of  the  forces  on  the  new  displacements  of  Eq.  41: 

Fli  =  Fli(ui,U3) 

Fl2  =  Fl2(u2,U3) 

FGi  *  Fgi(ui,U3) 

FG2  “  Fg2(u2,U3). 


(43) 


Fig.  10.  Displacement  components  of  three-degree-of- 
freedom  vehicle 


Ui  a  Uio  +  Vi 


i=l,2 


The  linear  approximation  of  lift  and  guidance  forces  can  be  expressed  as 
V  -p  I  .  dPLl.. 


F  -F  ■  ^^^2  ..  .  dFL2.. 
Fl2-Fl2o+-^V2+-^V3 


Fo,=FG,0+^n  +  ^V3 


Fg  =Fg  +— ^V2+^^V3 
'■02  "^^20  ^  ^  dV3  ® 

Using  Eqs.  39  and  45,  we  can  rewrite  Eq.  42  as 
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[M]{v}  +  [C]{v}  +  [K]{v}  =  0 


where 
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and 


C 

Cii  =  Ci2  =  — 


E 

C21  =  -C22  =  -  Y 


C33 


kii  =  - 


2  ^Fl^ 
m  avi 


(49) 


31 


ki2=- 


m  dv2 


ki3  =  “ 


m  ^  dv3  dv3 


,  2abdFGi 

I  dvi  I  dvi 


I  dv2  I  dv2 


k23=Mf!!Gi+!!^ 

“  I  I,  3v3  av3 


2b=' 


f!!k^!!k'l 

dv3  dv3  ^ 


k31  =  - 


k32=- 


1  ^FQi 
m  dvj 

1^FG2 

m  dv2 


k33  =  ~ 
m 


dv3  av3 


where 


9Fl. 

5Fg, 

^Fl,  ,  .  . 

3Fgo 

(49) 

(Cont’d) 


(50) 


32 


k//(h),  k/g(g),  kg/(h),  and  kgg(g)  are  motion-dependent  magnetic-force  coefficients 
(see  Fig.  5). 

If  we  assume  that  the  damping  effects  can  be  neglected,  the  eigenvalues  of 
Eq.  47  can  be  obtained  from 

[K]{v}=X[M]{v}  (51) 

where  X  =  cor  +  i  coi. 

With  magnetic  forces  and  stiffnesses  measured  by  the  experiments  (see 
Figs.  4  and  5),  the  eigenvalues  and  eigenvectors  of  a  maglev  vehicle  on  a  double 
L-shaped  guideway  were  calculated  with  the  theoretical  model  developed  in  this 
section.  Some  very  interesting  results  were  obtained  from  those  calculations. 

Figure  11  shows  that  eigenvalues  of  vehicle  motion  versus  levitation  height 
vary  when  guidance  gaps  are  fixed  (gj  =s  g2  =  Y*  =  12.7  mm).  The  first  mode  ©i 
shows  an  tmcoupled  heave  motion;  its  imaginary  part  of  the  eigenvalue  is  zero, 
while  the  second  and  third  modes  are  coupled  roll-sway  motions.  Within  the 
range  of  height  h  =  19.0  to  35  mm,  the  imaginary  parts  of  the  eigenvalues  appear 
not  to  be  zero.  This  indicates  that  within  this  range,  flutter  does  exist  for  these 
coupled  roll-sway  vibrations.  Table  1  and  Fig.  12  give  eigenvectors  and  modal 
shapes  of  these  three  modes  of  vehicle  motion,  respectively.  When  the  guidance 
gaps  are  fixed  at  gi  s  g2  Y*  s  5  mm,  the  same  results  are  obtained,  as  shown  in 
Fig.  13;  there  is  a  flutter  for  coupled  roll-sway  modes. 

Table  1.  Eigenvectors  of  vehicle  motion  (Y*  =  12.7  mm) 


h  = 

15.0  mm 

h 

=  25  mm 

h 

=  37  mm 

Mode 

VI 

V2 

V3 

VI 

V2 

V3 

VI 

V2 

V3 

Uncoupled 
heave  mode 

1 

1 

0 

1 

1 

0 

1 

1 

0 

©1 

Coupled 
roll-sway 
mode  ©2 

1 

-1 

-0.009 

0.586 

-0.586 

-0.332 

-1 

1 

-0.205 

Coupled 
roll-sway 
mode  ©3 

-0.545 

0.545 

1 

-0.810 

0.810 

0.060 

1 

-1 

0.448 

part  of  eigenvalue 
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Height,  mm 


Fig.  11.  MagleV'System  eigenvalues  vs.  vehicle  levitation  height, 
with  Y*  =  12.7  mm 


Uncoupled  Heave  Mode 


12.  Modal  shapes  of  three-degree-of-freedom 


I  mag 


Height,  mm 


Fi^.  13.  Maglev-system  eigenvalues  vs.  vehicle  levitation  height, 
with  Y*  -  5  mm 
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Figures  14  and  15  show  eigenvalues  of  vehicle  motion  versus  lateral  location 
of  vehicle  when  gi  =  g2  =  So  =  25  mm  and  levitation  heights  h  =  12.7  mm  and  h  = 
7  mm,  respectively.  We  notice  that  for  the  third  mode,  which  presents  the 
transversal  motion  of  vehicle,  the  real  part  is  zero  and  imaginary  peirt  is  not  zero 
within  a  certain  region.  This  indicates  that  the  divergence  is  subjected  to  the 
lateral  motion  of  the  vehicle,  given  those  vehicle  and  guideway  parameters. 
Figure  16  shows  the  real  part  of  the  third  mode  versus  lateral  location  of  vehicle 
when  the  parameter-equilibrium  guidance  gap  varies  as  gi  =  g2  =  go  =  10 
15  mm,  20  mm,  and  25  mm.  We  found  that  divergence  appesirs  only  in  the  case  of 
go  =  25  mm. 

We  must  point  out  that  the  measured  and  calculated  data  for  motion- 
dependent  magnetic-forces  and  force  coefficients  are  very  limited  and  that 
damping  effects  were  not  considered  in  the  above  analysis.  Even  though 
divergence  and  flutter  appear  in  the  eigenvalue  results,  we  still  have  difficulty  in 
completely  predicting  the  dynamic  instability  of  this  three-degree-of-freedom 
maglev  vehicle  model.  Further  research  is  needed  in  modeling  to  gain  an 
understanding  of  d3mamic  instability  in  maglev  systems. 


5  Closing  Remarks _ 

*  Motion-dependent  magnetic  forces  are  the  key  elements  in  modeling 
and  understanding  dynamic  instabilities  of  maglev  systems.  At  this 
time,  it  appears  that  very  limited  data  are  available  for  motion- 
dependent  magnetic  forces.  Efforts  will  be  made  to  compile 
analytical  results  and  experimental  data  for  motion-dependent 
magnetic  forces.  When  this  work  is  completed,  reconunendations 
will  be  presented  on  research  needs  on  magnetic  forces.  In  addition, 
specific  methods  to  obtain  motion-dependent  magnetic  forces  will  be 
described  in  detail. 

•  Various  options  can  be  used  to  stabilize  a  maglev  system:  passive 
electrodynamic  primary  suspension  damping,  active  electrod5mamic 
primary  suspension  damping,  passive  mechanical  secondary 
suspension,  and  active  mechanical  secondary  suspension.  With  a 
better  understanding  of  vehicle  stability  characteristics,  a  better 
control  law  can  be  adopted  to  ensure  a  high  level  of  ride  comfort  and 
safety. 


Imaginary  part  of  eigenvalue 


-30  -20  -10  0  10  20  30 

Lateral  location  of  vehicle,  mm 

Fig.  14.  Maglev -system  eigenvalues  vs.  lateral  location  of 
vehicle,  with  h  =  12.7  mm  and  go  =  25  mm 


Imaginary  part  of  eigenvalue  eigenvalue 
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Lateral  location  of  vehicle,  mm 


-30  -20  -10  0  10  20  30 

Lateral  location  of  vehicle,  mm 

Fig.  15.  Maglev -system  eigenvalues  vs.  lateral  location  of 
vehicle,  with  h  =  7  mm  and  gg  =  25  mm 


Lateral  location  of  vehicle,  mm 


Fig.  16.  Real  part  of  maglevsystem  eigenvalues  vs.  lateral 
location  of  vehicle,  with  h  =  7  mm  and  go  =  10,  15, 

20,  and  25  mm 

*  Computer  programs  are  needed  to  screen  new  system  concepts, 
evaluate  various  designs,  and  predict  vehicle  response.  It  appears 
that  the  stability  characteristics  of  maglev  vehicles  under  different 
conditions  have  not  been  studied  in  detail  in  existing  computer  codes. 
When  information  on  motion-dependent  magnetic  forces  becomes 
available,  the  existing  computer  codes  can  be  significantly  improved. 

*  Instabilities  of  maglev-system  models  have  been  observed  at  Argonne 
and  other  organizations.  An  integrated  experimental/anal3rtical 
study  of  stability  characteristics  is  an  important  aspect  of  maglev 
research. 
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Fig.  9.  Maglev  system  with  a  vehicle  operating  on  double  L- 
shaped  aluminum  sheet  guideway 

a  =  0.5  H  (nun) 

b  =  0.5(W  -  25.4)  (mm). 

Equations  of  motion  for  this  three-degree-of-freedom  maglev  system  can  be 
written  as 

m2  +  Cz=  Fl^  +Fl2  -  mg 

I0  +  E0  =  a^FQj +FQ2)  +  b(FL^ -Fl^)  (38) 


my  +  Dy  =  FG^+FG2, 


where  m  is  the  mass  of  the  vehicle,  C  and  D  are  damping  ratios;  I  is  the  moment 
of  inertia  about  the  center  of  mass  of  the  vehicle  [I  =  (m/12)(H2  +  W2)].  Fl^,  Fl2. 
Fgj,  and  Fg2  sire  lift  and  guidance  forces  and  are  functions  of  y  and  z.  At  the 
equilibrium  position,  they  are  FLio(yo»Zo)»  FL2o(yo.Zo),  FGio(yo,Zo),  and  FG2o(yo.Zo)- 
Apply  them  to  Eqs.  38: 
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FLio  =  FL20 
Flio  +  Fl20  =  nig 
FQio  “  ~Fg2o 
Therefore 

m  ~  ^Lio  "^FLao  _  2FL(hQ) 

g  g 

Let 


(39) 


(40) 


z  =  |(ui  +  U2) 

y  =  ^3  (41) 

9  =  (ui  -U2)/2b 

j 

where  u^,  U2,  and  113  are  shown  in  Fig.  10.  Equations  38  can  be  rewritten  as 
m(ui  +  u2)  +  C(ui  +  u2)  =  2(Fl^  +Fl2  ~nig) 

l(ui  -  U2)+  |(ui  -  U2)  =  2a(FG,  +  Fg2  )  +  2b(Fi,,  -  Fl^  )  (42) 

mu3  4Du3  =  FGj+FG2. 

Note  the  reduced  dependence  of  the  forces  on  the  new  displacements  of  Eq.  41; 

Fli  =  Fli(ui,U3) 

Fl2  =  Fl2(u2,U3) 

FQi  =  Fgi(ui,U3) 

Fg2  =  Fg2(u2,U3). 


(43) 
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Fig.  10.  Displacement  components  of  three-degree-of- 
freedom  vehicle 


Ui  =  Uio  +  Vi 


i=l,2,3 


The  linear  approximation  of  lift  and  guidance  forces  can  be  expressed  as 

Fg  =Fg  +^55Xvi  +  ^^V3 
^Gi  ^Gio  1  3 

rn  „  3Fg2 

Fg3=Fg,„+-^V2+-^V3 


Using  Eqs.  39  and  45,  we  can  rewrite  Eq.  42  as 
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..  ..  C.  C.  2  3Fli  2  3Fl2  2f3FL,  3Fl  \ 

V1  +  V2+  — Vi  +  — V2 - - r-^V2 - +  V3=0 

m  m  m  dvi  m  dV2  dv3  dv3  I 


..  ..  E  .  E  .  f2ab  3Fg,  2h^  ^Fl,  ' 

■Vl  +  V2--v,^-V2+  — ^  ^  V, 


^  r  2ab  3FGg  2b^  ^Fl^ 
I  I  3v2  I  dv2 


I  [  dv3  dv3  1  I  [  dv3 


V3  =0 


.D.  1  ^Fg^  1  ^Fg2  if^FG^  9Fg2\ 

V3+— V3 - - r— ^V2 - +  V3=0 

m  m  dvj  m  dv2  ml^  dv3  dv^  i 


[M]{v}  +  [C]{v}  +  [K]{v}  =  0 
where 


■  1  1  O' 

Cii  C22  0 

'kii  ki2  kia' 

M  = 

-1  1  0 

,  c  = 

C21  C22  0  |,  K  = 

^21  ^22  1^23 

0  0  1 

.  0  0  C23J 

.^31  ^32  k33 

Cll  =  Cl2  =  — 
m 


C2I  =  ~C22  =  “Y 


D 

C33  =  — 
m 


m  dvj 


(49) 


m  3v3  3v3  J 

2ab  ^  2b^ 

I  9vi  I  3vi 

2ab  2b^  3Fl2 

I  ^V2  I  dv2 


2ab 

r3B^ 

^aFo,] 

2b^ 

I 

1  ^V3 

dv3  J 

I 

^  ^V3 

av3  j 

1 

m  3vi 
1 

m  dv2 


(49) 

(Cont’d) 
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